RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 69D, No. 2, February 1965 

Momentum Transfer Collisions in Oxygen 
for Thermal Electrons 

Michael H. Mentzoni 

Contribution From the Applied Research Laboratory, Sylvania Electronic Systems, a Division of Sylvania 

Electric Products, Inc., Waltham, Mass. 02154 

(Received August 10, 1964; revised September 25, 1964) 

The electron collision frequency for momentum transfer, v m , has been determined by 
microwave methods in an oxygen plasma for temperatures between 300 and 900 °K. It 
was found that the measured average electron-neutral collision frequency was consistent 
with vem being proportional to the square of the electron velocity when a Maxwellian velocity 
distribution is assumed for the electrons; e.g., „ em ( w ) = 4.4X 10~ 8 X iV(0 2 ) x u (sec" 1 ), where 
u is electron energy in electron volts and N (0 2 ) is the number density of oxygen molecules. 
The results also indicated that electron-ion collisions are important at the Lower temperatures 
with the values of v ti in agreement with the ones recently reported by Chen. 



1. Introduction 

The electron collision frequency for momentum 
transfer, v m and its energy dependence in atmos- 
pheric gases is of special technical interest. 

In atmospheric environment the total electron 
collision frequency for momentum transfer is deter- 
mined by adding the collision frequencies for the 
various gas constituents in proportion to their 
number densities. Thus by laboratory determina- 
tions of v m as a function of electron energy for gases 
such as N 2 , 2 , NO, N, O, one is able to compute 
the composite value of v m . This value will vary 
with altitude according to the density profiles of the 
species and the electron temperature, the latter 
covering the temperature range of the present investi- 
gation between altitudes from about 110 to 150 km 
[Nawrocki and Papa, 1961]. 

The collisional absorption and its effects on earth- 
space radio propagation has been recently reviewed 
by Lawrence et al. [1964] who, in this respect, 
points out the importance of D layer, the physical 
processes of which have been surveyed by Reid [1964]. 
Calculations of EM-wave propagation based upon 
the Appleton-Hartree magneto-ionic theory incor- 
porate an effective velocity independent collision 
frequency thus simplifying the analysis. It has 
been shown by Phelps [i960], that such a simplifica- 
tion may lead to significant errors. 

During reentry, signal absorption due to the 
ionized shock layer and wake constitutes a serious 
problem the solution of which necessitates a knowl- 
edge of v m in shock heated air [Bachynski et al., I960]. 



Temperatures in a wake are in general less than 
3,000 °Kfor reentry speeds and the v^ measured here 
would pertain to the far wake which is of importance 
for radar detectability 

A number of laboratory studies have been reported 
in recent years pertaining to nitrogen. For oxygen 
and dry air the available information is still rather 
limited [Nielsen and Bradbury, 1937; Healey and 
Reed, 1941; Van Lint et al./ 1959; Phelps, 1960; 
Shkarofsky et al., 1961; Carruthers, 1962; Gilardini, 
1963]. For thermal and near thermal energies, the 
dependence of v m on the electron energy u has been 
difficult to determine experimentally. This is so 
because the generation of monoenergetic electron 
beams, used so successfully in scattering experi- 
ments above 1 eV becomes increasingly difficult at 
lower energies. However, Pack and Phelps [1959] 
have improved the design of the mobility tube 
allowing measurements at thermal energies. It 
seems that the reanalysis by Shkarofsky et al. [1961] 
of results from the early drift tube experiments of 
Nielsen and Bradbury [1937] with oxygen pertaining 
to u higher than the gas temperature, T gaa , is con- 
sistent with a collision frequency given by the formula 

v em N(0 2 )- 1 =^XlO- 9 u 1/2 +5XlO- 8 u cm 3 sec" 1 



where iV(0 2 ) 
oxygen and 



is the molecular number denisty of 
m refers to elastic collisions between 
electrons and molecules. Measurements pertaining 
to T ga8 = 300 °K, assuming thermalized electrons, 
give results that are in agreement with the v em N(0 2 )~ l 
= 7Xl0~ s u cm 3 sec" 1 (Phelps, 1960); it seems 
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desirable to perform direct measurements at some 
other gas temperatures within the suggested range 
of applicability of this equation and thus obtain a 
more complete verification. 

Finally, electron-ion collisions, v ei , which are 
important at high electron-ion densities and lower 
temperatures are considered and their collision cross 
section evaluated. For a review of v em and v ei per- 
taining to the atmosphere the reader is referred to 
publications by Nicolet [1959], Nawrocki and Papa 
[1961], and Dalgarno [1961]. 



2. Experimental Method and Procedure 

The experimental arrangement has been described 
previously [Mentzoni, 1963; Mentzoni and Row, 
1963] and only the salient points will be discussed 
here. 

A short duration, pulsating d-c voltage, periodi- 
cally created an oxygen plasma inside a cylindrical 
quartz tube which also functioned as a waveguide 
operating in the "quasi" TE n mode due to a gold 
coating on the outer surface. A low power X-band 
microwave signal propagating along the 80 cm long- 
plasma column acted as the probing wave. An 
interferometer technique together with accurate 
measurements of plasma insertion losses permitted 
determination of the ratio. 
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here AA, A(j>, Aa, and Aft represents changes in the 
total attenuation (nepers), phase shifts (radians), 
attenuation and phase constant respectively due to 
the plasma, L is length of the interaction region, 
E is electric field intensity, R is inner radius of the 
quartz cylinder and <r r and o- t are the real and imagi- 
nary components of the rf conductivity 
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Here N e , e, m, and v are electron number density, 
charge, mass, and velocity, respectively, and f is 
the spherically symmetrical part of the electron 
velocity distribution function normalized such that 
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Equation (1) is valid provided the first order pertur- 
bation analysis is sufficiently accurate [Golant and 
Zhilinskii, 1960] which means that 
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A F goes from unity to 1.3 for the TEn mode as N e 
goes from a uniform to a J (2A05r/R) distribution. 
Actually the quartz walls cause a further (and in 
this case insignificant) perturbation of the TE n 
mode as shown by Unger [1957]. 

Even though the probing field is nonuniform, it 
is selected to be very weak, satisfying the criterion 
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with (u) and G the average electron energy and 
energy loss factor (G— factor), respectively, so that 
one has that the average electron energy is much 
larger than the energy perturbation caused by the 
probing field, i.e., 
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and thus can be regarded as spatially uniform. 
When a larger field intensity is used, such as in the 
cross modulation experiment, the distribution func- 
tion can still be regarded as independent of position 
provided the thermal conductivity of the electrons 
is very large. 

In the present arrangement the gas is heated to 
any desired ambient temperature up to about 1000 
°K by surrounding the waveguide-discharge con- 
figuration by an electric oven. A microwave radiom- 
eter utilizing the noise sampling technique is used to 
monitor the electron radiation temperature, one 
can thus find the time after which isothermal condi- 
tions prevail during the afterglow [Mentzoni, 1964] 
which ranged typically from about 30 to 50 jusec. 
In comparison, the total measurable time decay for 
electrons and ions, characterized by iVe>5X 10 8 cm" 3 , 
lasted always more than 1 msec. 

In general one has v m =v em -\-v ei J rv ee , where the 
new right-side terms indicate the frequencies for 
electron-ion and electron-electron collisions. In 
the following v ee will not be considered. Whereas 
v em is independent of spatial coordinates, this is, in 
general, not so with v ei which depends on the ion 
density, N t according to the relation [Shkarofsky, 
1961]. 

v ei =C 1 ^hi[C 2 Tr 2 N7^] 

where C x and C 2 are constants. 

Ignoring the logarithmic term one can write, 
since Ni=N e , 

Vei=-V°ei(v)F(r) 
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where F(r) describes the electron density profile 
defined by N e (r)=N F(r), with N being the electron 
density at the center. 

When isothermal conditions prevail (T e =T g r lB ) 

the electrons obey a Maxwellian distribution so 

3 
that (u)=-kT gas , this is also a reasonable assump- 

tion when T e is slightly in excess of T gas due 
to an rf heating pulse (cross modulation). In 
this case and for v m <C<Cu (low pressure approxi- 
mation) one ha; from eqs (1) and (2) 
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It is seen that for a nearly uniform electron distribu- 
tion, F(r)~l one has 
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In general v m is a function of electron energy u(= 
jimv 2 ) and it can be shown [Phelps et al., 1951] 
that one can deduce from the measured average 
quantities the energy dependence of v em under 
certain conditions. In this study the condition 
v 2 m <^<iu) 2 is satisfied and the right side of (3) can 
be written, assuming v em ^>^>v e i, 
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Based on previously reported results it is reasonable 
to assume a power law v em =NCv h and one obtains 
from (4) 
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where G and h are constants which may be found 
empirically from a best fit procedure, and N is 
gas density in molecules per cubic centimeter. 
If indeed there exists such a simple power law for 
v em and measurements are performed with the same 
N, but various temperatures, T e =T gAa , one obtains 
h from the simple relation 



(Vem)Tek = (&A/A(l>) Tek _ 
(v em) Tel (AA/A<t>) T el 



m 
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If h can be determined in this manner, G can be 
determined from (3) and (4). 
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Figure 1. — Typical lime variations of <Cv m ^> as computed 
from plasma insertion-loss measurements during the after- 
glow. 

The scatter of the experimental points represents differences between several 
runs of measurements. The broken lines Indicate where the steady-state value 
of Om> should bave been, based on the one found at /;=3.0 torr. 



3. Results and Discussion 

In the preceding section it was shown how it is 
possible to obtain v em {u) for the low pressure approxi- 
mation under the assumption that v em depends on v 
according to a simple power law. Figure 1 shows a 
plot of (A yl/A0)co versus time at r gaB =894 °K. It 
is seen that isothermal conditions are not present 
during the very early afterglow, thus perhaps 
invalidating (4). This is so because the assumption 
J being Maxwellian may not be true even though 
the inequality j> m «« (=2t 10 10 cycles sec _1 ) is 
well satisfied. At later times constancy of (A A/ 
A0)o? is obtained indicating T e =T gQi% in agreement 
with electron radiation temperature measurements. 
The results for this case are given in table 1 and it 
is seen that h=2 yields good agreement with experi- 
mental results. The relationship thus found reads 
v em (u)=4:.4X lO~ 8 Xu N{0 2 ) sec -1 with u in electron 
volts or <> cm > = 3.0Xl0 5 pT e sec" 1 , with p (torr) 
referred to 273 °K. At temperatures below 569 °K 
steady state values of Om> as found for the higher 
temperatures were not as conclusive pointing to 
the possible noticeable contribution from <*> 6 *>. 

The spatial distribution of Ni and thereby v ei 
must in this connection be considered since (3) is 
the basis for the interpretation of the measurements. 
As has been pointed out [Mentzoni, 1965] the 
plasma seems to be axially and tangentially uniform 
based on luminosity measurements. The radial 
variation of N e (r) depends on the initial distri- 
bution and the rate of ambipolar diffusion during 
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the afterglow. At the pressures pertaining to the 

present experiment, it is reasonable to assume a 

radially uniform initial distribution. The ratio 

C R C R 

E(r) 2 F(r) 2 r dr/\ E(r) 2 F(r) r dr appearing in (3) is 
Jo Jo 

the spatial average F(r) of the electron-ion profile 
weighted by the quantity E{r) 2 F(r) with E(r) vary- 
ing approximately as jAk-^ ), thus, F(r) is most 
sensitive to F{r) in the center region. 

Table 1. — Ratios of ^^m^ for various temperatures under 
isothermal conditions compared with the case of *> em = NC 
v b (h = 2) and <> em > being averaged over a Maxwellian 
velocity distribution 



Tek 


Tel 


Kvem^Tek/Kvem^Tel 


(Tek/Tel) h/2 

h=2 


894 
734 
894 
569 


734 
569 
569 
450 


1.27 

1.29 
1.64 
1.4 


1.218 
1. 290 
1.571 
1.264 



Table 2. Electron-ion density profiles for various times during 
the afterglow, T = 300 °K, p = 0.3 torr and 1)^ = 110 cm 2 
sec* 1 torr l 
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10 


0.980 


0.980 


0.980 


0.980 


0.961 


20 


.980 


.980 


.980 


.979 


.880 


30 


.980 


.980 


.980 


.973 


.796 


40 


.980 


.980 


.980 


.959 


.725 


50 


.980 


.980 


.978 


.959 


.667 


60 


.980 


.980 


.975 


.915 


.618 



1 D a is am bipolar diffusion coefficient pertaining to the ions O^. 



Table 3. Some composite results of <^ m ^>= r <^em>4-<C» / ei> 
for T gaa = T e = 300 °K 
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A = \0-*X<v tm >° 


B=10-*X<l>ei> b 


A+B 


10-SX<^> C 


cm-3 


Torr 


sec- 1 


sec- 1 


see- 1 


sec- 1 


3.325 
6.650 
9.975 


3.0 
3.0 
3.0 


2.70 
2.70 
2.70 


1.53 
2.89 
4.20 


4.23 
5.59 
6.90 


4. 35±. 85 
6. 07±L 17 

8.5 


6.650 
9.975 


1.5 
1.5 


1.35 
1.35 


2.89 
4.20 


4.25 
5.55 


3. 65±. 05 
5.0 



° Based on present results, <t>em>=3.0XpXT e X10 5 (sec- 1 ). 



» Based on Chen's experimental results, O et >=3.6 
(sec -1 ). 
c Present direct measurements. 
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t Considering diffusion losses under isothermal con- 
ditions only [Mentzoni, 1964] F(r) can be computed 
as illustrated in table 2. At room temperature, 
y g as=300 °K, some of the results obtained are illus- 
trated in table 3. Due to the lack of isothermal 
conditions in the time range of the high electron-ion 



densities at the lowest pressures the experimental 
points tend to lie below those to be expected from 
theory and previous experiments (see fig. 2) assuming 
r e =T gas =300 °K which in general follow the law 



<^> = C^T-^m {C 2 TTN-\> 2 



(6) 



where Oi, C 2 are constants. Several theoretical and 
two experimental evaluations of 0\ and C 2 have been 
made previously yielding values of d which varies 
between the limits 2.25 and 4.84X10 3 , and C 2 
between 3.32 X10 3 and 20X101 

It should be noted that Ginzburg and Gurevich 
[1960] and Nicolet [1959] quote expressions for 
<C*^v> which are somewhat modified in form as 
compared to (6) being <i> 6f >~5.5iV*2V 3/2 In 
(2Xl0 2 T e /N e v*) and <^> = [59.1+4.18 log 10 (T// 
N e )] Ni/T e V2 10~ 6 , respectively, in the latter equation 
the electron density is given per meter. 3 At the pres- 
sures ^=3.0 mm Hg and p=1.5 mm Hg it is seen 
from table 2 that the agreement with Chen's [1964] 
results is fair considering the experimental error of 
the measurements. 

From figure 2 it is seen that the experimental 
points for the two different pressures converge for 
the lower electron densities and electron temperature 
indicating that <^ e *> is dominating. At the higher 
densities on the other hand, the increasing divergence 
of the experimental points for the two pressures is 
to be expected in the case <Cy em ^> is dominating. 
The solid curve in figure 2a indicate that even under 
nonisothermal conditions the present result for 
<C^™> coupled with <Cv ei ^> as found by Chen is 
reasonably consistent with experimental observation 
of <^v m ^> and electron radiation temperature meas- 
urements. As the electron energy increases the 
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Figure 2. — (b) <v m > versus the differential phase-shift 
A<j> for T eaa = 300 °K and pressures 0.30 and 0.15 torr. 

The measured values are indicated by points, the dash-dotted curve repre- 
sents <i>ei> under isothermal conditions according to Chen, and the solid curve 
is computed from <v m >=<rem>+<j>,i>=S.O p T e X10*+3.6 
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with p = 0.30 torr (referred to 273 °K) and T e following the time decay curve 
shown in (a) . The solid portion of this curve (time>20 /xsec) has been measured. 
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inelastic electron collisions involving rotational and 
vibrational excitation of the oxygen molecules will 
also become increasingly important, however, their 
contribution to <Cv m ^> is expected to be small at the 
temperatures considered here. 
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